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Abstract 

The connection between ghost-free formulations of RG-invariant perturbation theory in the 
both Euclidean and Minkowskian regions is investigated. Our basic tool is the "double spectral 
representation", similar to the definition of Adler function, that stems from first principles of 
local QFT. It relates real functions defined in the Euclidean and Minkowskian regions. 

On this base we establish a simple relation between 

— The trick of resummation of the 7r 2 -terms (known from early 80s) for the invariant QCD 
coupling and observables in the time-like region and 

- Invariant Analytic Approach (devised a few years ago) with the "analyticized" coupling 
d&n{Q 2 ) and nonpower perturbative expansion for observables in the space-like domain which 
are free of unphysical singularities . 

As a result, we formulate a self-consistent scheme — Analytic Perturbation Theory (APT) 
- that relates a renorm-invariant, effective coupling functions a an (Q 2 ) and a(s) , as well as 
non-power perturbation expansions for observables in both space- and time-like domains, that 
are free of extra singularities and obey better convergence in the infrared region. 

Then we consider the issue of the heavy quark thresholds and devise a global APT scheme for 
the data analysis in the whole accessible space-like and time-like domain with various numbers 
of active quarks. 

Preliminary estimates indicate that this global scheme produces results a bit different, some- 
times even in the five-flavour region, on a few per cent level for a s - from the usual one, thus 
influencing the total picture of the QCD parameter correlation. 
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1 Introduction 



1.1 Preamble: perturbation theory and a s 

The issue of the strong interaction behavior at the low and medium energy W = y/s and 
momentum transfer Q = y/Q 2 attracts more and more interest along with the further 
experimental data accumulation. 

As a dominant means of theoretical analysis, here one uses the perturbative QCD 
(pQCD), in spite of the fact that in the given domain the power expansion parameter a s 
is not a "small enough" quantity. Physically, this region corresponds to three (/ = 3) 
and four (/ = 4) flavors (active quarks). Just in the three-flavour region there lie 
unphysical singularities of central theoretical object - - invariant effective coupling d s 
. These singularities, associated with the scale parameter Aj =3 ~ 350 MeV, complicate 
theoretical interpretation of data in the "small energy" and "small momentum transfer" 
regions Q = \J~Q 2 ~ < 3A 3 .) On the other hand, their existence contradicts some 

general statements of the local QFT. 

It is important to notice that in the current literature for the effective QCD coupling 
in the time-like domain a s (s) ; s = W 2 one uses literally the same expression, like one 
in the Euclidean domain. By the way, implanting the mentioned singularities into the 
three-flavor region of small energies W ~ 350 MeV. 

Meanwhile, the notion of invariant electron charge (squared) a{Q 2 ) = e 2 (Q) in QED 
has initially been defined in the early papers|J on renormalization group (RG) only in the 
space-like region in terms of a product of real constants Zi of finite Dyson renormalization 
transformation. Just the Euclidean invariant charge e(Q) is related by the Fourier trans- 
formation to the space distribution e(r) of the electric charge (around a point "bare" 
electron) introduced by Dirac0. 

Analogous motivation in the RG formalism (for detail, see chapter "Renormalization 
group" in the text||) underlies a more general notion of invariant coupling g(Q) , defined 
only in the space-like domain. Inside the RG formalism, there is no simple means for 
defining g in the time-like region. 

Nevertheless, in modern practice, inspired by "highly authoritative reviews" 0, |J and 
some monographs (like ||) one uses the same singular expression for the QCD effective 
coupling a s both in the space- and time-like domains. 

Technically, this "implanting" of the a s Euclidean functional form into Minkowskian 
is accompanied by some modification of numerical expansion coefficients. To the initial 
coefficient calculated by Feynman diagrams, one adds specific terms (containing it 2 and its 
powers) with coefficients of some lower orders. These 'V 2 -terms" are the only "atonement 
for the Styx river crossing" from the Euclid realm to the Minkowski domain. 

1.2 Time— like region, it 2 — terms 

Meanwhile, as it easy to show, the "7r 2 -procedure" is valid only at small parameter 
7r 2 / ln 2 (s/A 2 ) values, that is in the region of high enough energies W ^> Ae 71 "/ 2 ~ 3 GeV . 
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Here, it is useful to restore the ideas proposed by Radyushkm0 and Krasnikov — 
Pivovarov|| (RKP procedure) at the beginning of the 80s. 

To introduce an invariant a s in the time-like region, both the authors used integral 
transformation R, "reverse" to the Adler function definition. The last one can be treated 
as the definition of integral operation 



R(s)^D(z) = Q 2 T ^-R(s)^B{R(s)} 
Jo V s i z ) 



1) 



transforming a real function R(s) of positive (time-like) argument into the function D(z) 
defined in the cut complex plane with analytic properties adequate to the Kallen-Lehmann 
representation. In particular, D(Q 2 ) is real at the positive semiaxis z = Q 2 + i0 ; Q 2 > . 
The reverse operation R can be expressed via the contour integral 



R(s) 



2tt 



i+ie 



dz 



D pt (-z) = R [D(Q' 



By operation R one can define RG-invariant, effective coupling a(s) = R [a s ((5 2 )] in 
the time-like region. A few simple examples are in order 

- For the one-loop effective coupling ai 1 -* = [f3 ln(Q 2 /A 2 )] -1 one has || 0, p~0|| F C l 



Po 



1 1 L 
arctan — 

2 7T 7T 



L>0 



1 7T 

— — arctan — 

POTT L 



At the two-loop case, to the popular approximation 



A)«S op (Q 2 ) = y 



h(f) /2 



In/ 

J 2 



I = In 



A 2 



there corresponds J7|, |TTJ 



L = ln-. (2) 



a pop( S ) 



biL 



L 2 + 7P 



bt In [V^ 



7T Z 



+ 1 



A) 



L 2 + 7T 2 



(3) 



Both the expressions (|2]) and (0) are monotonously decreasing with finite IR limit 
fi(0) = l/A)(/ = 3)~1.4. 

- At the same time, square and cube of ats transform into simple "pipizated" 
expressions^, § 



a«( s ) = r 



which are not powers of a^\s) . 



P 2 [L 2 + vr 2 



(31 [L 2 + ir 



i2 ' 



(4) 



1 Later on, this idea has been discussed by several known authors — see Refs. [12 ■ [fL5[. 
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Note also that transition from singular a s and its powers to "pipizated" expressions, 
that is operation R, can be performed [T^ , by the differential operator 

sin 7r7^ d 

TZ = — — — ; V = Q 2 — — tt with the substitution Q 2 — > s , e.g., 7Za s (Q 2 ) = a(s) . 



The most remarkable feature of all presented expressions for a(s) and Slfc(s) (valid in a 
more general case) is the absence of unphysical singularity (the "log pole" at the one-loop 
case)0 which is "screened" by tt 2 -contributions. 

Besides, a common "Euclidean" perturbation expansion 

D pt (Q 2 ) = l + J2d k a k s (Q 2 ) (5) 

k>l 

in powers of the standard RG-summed effective coupling a s (Q 2 ), with its unphysical 
singularities in the IR region (at Q 2 < A 2 ,) being transformed by R to the time-like 
region, transits into the asymptotic expansion over a nonpower set of functions 

R^s) = R [D pt (Q 2 )] = 1 + 52d k K k (s) ; = R [a k s (Q 2 )] , (6) 

k>l 

with better properties of decreasing |7[ of subsequent terms . 

At the same time, higher functions, like (^), vanish 21^(0) = ; k > 1 in the IR limit. 

On the other hand, in the UV region at ln(s/A 2 ) ^> n , i.e., for W 3> Ae 71 "/ 2 ~ 
3 GeV , the functions a and can be represented as a series in powers of the parameter 
7r 2 /L 2 , L = ln(s/A 2 ) . Such expressions sometimes can be reformulated into expansions 
in powers of a s . For instance, in the one-loop case 

- jj; ~ ^ + ^+ - ~ W\s)) 3 + ^ (a?\s)) 5 (7) 



Without going into detail, note?] that, qualitatively, the functions 21^ behave very 
similarly to the functions A k involved into non-power Euclidean asymptotic expansion for 
observables||18|| arising in the Analytic Perturbation Theory (APT) - see below Section 
p..3| and Figure 2. In particular, they oscillate at small argument values and form an 
asymptotic set a la Erdelyi. 

As it follows from eqs. (||) and @), one-loop "pipizated" functions satisfy the recur- 
sion relation (d / dL)W k {s) = — k /?o 2lfe+i( s ) which is analogous to the one-loop differ- 
ential equation for the invariant coupling. According to |nj, this recursion is valid for 
analyticized functions A k ^. 



2 This feature was not mentioned in the pioneer papers of the 80s we have cited above. 

3 See, also the first version of this paper ||16| . A more minute numerical information on functions a 



c*an , 212,3 and -42,3 can be found in recent paper by Magradze flTi 
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Quite recently |l7j the two-loop generalization has been found 

llL^ ){s) = " " ^U s ) ■ (8) 

It can be considered as (at k — 1 ) a mould of two-loop differential equation for a s . 
Analogous relation is valid for analyticized Ak ■ 

For the reverse transition from Minkowski to Euclid, one could try to use the trans- 
formation D defined by ([!]). However, it is evident that we shall not return to the initial 
coupling a s and to series in its powers (|). To elucidate the issue, it is useful to turn to 
the foundation of the Invariant Analytic Approach mentioned above. 



1.3 Space— like region : APT 

Indeed, as it has been well-known from the late 50s |L9| , there exists a method of getting rid 
of Euclidean unphysical singularities by combining RG-summed expressions with Kallen- 
Lehmann analytical representation for 6t s (Q 2 ) in the Q 2 variable. In the mid-90s this 



idea was used in QCD |2(J, ^T], |22[ under the name of Invariant Analytic Approach. Its 
further development and application to perturbative expansion for observables yielded 
Analytic Perturbation Theory — [f23] . 



We remind here the basic features and results of APT ( — see also a recent review 

til)- 

By combining three elements 

1. Usual Feynman perturbation theory for effective coupling(s) and observables, 

2. Renormalizability, i.e., renormalization-group (RG) invariance, and 

3. General principles of local QFT — like causality, unitarity, Poincare invariance 
and spectrality — in the form of spectral representations of the Kallen-Lehmann and 
Jost-Lehmann-Dyson type 

it turns out to be possible to formulate an Invariant Analytic Approach (IAA) for 
the pQCD invariant coupling and observables in which the central theoretical object is a 
spectral density. 



• Being calculated by the usual RG-improved perturbation theory, it defines and 
relates Q 2 -analytic, RG-invariant expressions for effective RG-invariant coupling 
and perturbative observables in the Euclidean channel. 



In particular, the IAA results in the modified ghost-free expression for the invariant 
QCD coupling a an (<5 2 ; /) which is free of ghost troubles and obey reduced [f2T 
higher-loops and renormalization-scheme sensitivity^. See, Fig.l. 



4 This analyticized QCD coupling a an has been successively used|29[ in analysis of the pion and 
7*7 — > 7r° formfactors. 
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The IAA change the structure of perturbation expansion for observables: Instead 
of common power series, as a result of integral transformation, there appears non- 
power asymptotic expansion |18 a la Erdelyi over the sets of specific functions 



Ak(Q 2 ', f) , free of unphysical ghosts. These functions are defined via integral trans- 
formations of related powers a^(Q 2 ; f) in terms of relevant spectral densities. This 
nonpower expansion for an observable, with the coefficients extracted from the rel- 
evant Feynman diagrams, we call the Analytic Perturbation Theory. 

At small and moderate arguments, At diminish with the k growth much quicker 
than the powers of a^ n (and even oscillate in the region -y/i, Q ~ A ,) thus improving 
essentially the convergence of perturbation expansion for observables. 

The first purpose of this work is to elucidate relation between the Radyushkin- 
Krasnikov-Pivovarov procedure leading to effective summation of 7r 2 -terms ( "pipization" 
trick) [^, Q for observables and the Solovtsov[|II], |T3| construction of the effective QCD 
coupling within the IAA scheme in the s-channel. 

In the course of this analysis — see Section [2] — we discuss the APT proliferation 
to the time-like region, remind a spectacular effect of "distorting mirror" correlation ||27|| 
between analyticized and pipizated invariant QCD couplings in space-like a an (Q 2 ; /) and 
time-like a(s; /) regions (see Fig.l below), and establish this effect for corresponding 
expansion functions Ak(Q 2 ', f) and 2lfc(s; /) - see Fig. 2. 

Then, in Section [3], we consider an the transition across the heavy quark thresholds, 
to construct a "global" picture for the whole physical region M T < Q < Mz — see 
Fig.2. 



It should be noted, that all precedent papers Refs. || - |36| dealt only with the massless 
quarks in the case with fixed flavour number / . This can be justified, to some extent, 
when analyzing inside a narrow interval of the relevant energy y/s or momentum transfer 
Q values. Meanwhile, the ultimate goal of all the pQCD is a correlation of effective 
coupling values extracted from different experiments. 

Main results of this investigation are reviewed in the Conclusion. 

2 Self-consistent scheme for observables 
2.1 Modification of the APT 

As it has been mentioned above, applying operation D to a does not restore a usual 
effective coupling as far as representation ([I]) is not compatible with ghost singularity of 
a s . 

Instead, we arrive at 

D {a(s; /)} = / a * a(s; f) = « an (Q 2 ; /) , (9) 

71 Jo V s + v ) 
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i.e., to effective Euclidean coupling a an of APT. This simple fact has first been established 
in [HI. We see that operations D and R relate "pipizated" and "analyticized" coupling 



functions in space- and time-like regions. Hence, in this case R = D -1 . Note, however, 
that the relation DR = 1 is valid only for the class of functions F(Q 2 ) e Ckl satisfying 
the Kallen-Lehmann representation. 

Now, we have the possibility of extending the APT to the time-like region. We shall 
do it in the form of a recipe, using operation of analyticization 

F(Q 2 ) -> F an (Q 2 ) = A ■ F(Q 2 ) , (10) 

first introduced in Refs. pO, |21| in terms of the Kallen-Lehmann representation 



oo 

FUQ 2 ) = ~ f-^n* = $*F(-a), (11) 

7T J 0- + Q Z 


with spectral density defined via straightforward continuation of F on the cut. 

Relations fllOD and (|TT|) together define the analyticization operation. Now, we can 
formulate the APT anew. 

Firstly, one has to transform the common singular coupling function a s (Q 2 ) or some 
power expansion of an observable 

D pt (Q 2 ) = l + J2d k a k s (Q 2 ;f) (12) 



fc>i 



into the corresponding analytic Euclidean expression a an or D an (Q ) , free of ghosts 

DUQ 2 ; /) = i + Yl dk Ak (Q 2 '> /) ' a »(^: /) = MQ 2 \ f) , (is) 



k>l 



oo 

A k (x; f) = - [ Pk (a; f) ; p k (a; f)=Q [a k s (-a; /)] 

7T J a + x 



(14) 



with spectral densities p , p k introduced according to (|IT 



Secondly, by operation R one defines fpTofl in the Minkowskian region invariant coupling 
function^ 



« an (Q 2 ; /) - 5(s; /) = R = f ^-p(a; /) (15) 

s 

or some other quantity like 

R n (s) = R [D pt (Q 2 )} = l+J2 d ^{s) ; Sl fc (s) = R K"(^ 2 )] ( 6 ) 



fe>i 



5 As it follows from this expression, the spectral function can be considered as a beta-function. How- 
ever, contrary to Schwinger's hope|37j], this p(s; f) , being a spectral function for the Euclidean invariant 
coupling, happens to be the RG generator for another, Minkowskian, invariant coupling |lq]. 
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with 



oo 



(16) 



Finally, we have a simple possibility of reconstructing an Euclidean object from the 
corresponding Minkowskian one with the help of the dipole operator D like 

a an (Q 2 ;/) = D {a(s; f)} . 

In particular, by substituting /) into the integrand, we obtain after integration by 

parts 

Q 2 f°° da fl 1 ln(o-/A 2v 



D{qW} 



o Jo 



1 



(a + Q 



+ 



2^2 



arctan 



TT 



7T 



A 2 



(17) 



_lnQ 2 /A 2 A 2 -Q 2 

This simple calculation elucidates the relation between ghost-free expressions in the 
Minkowskian and Euclidean regions. They are related by a reverse transformation as well. 
For instance, in accordance with (Il5|), 



5«i 



s;f) = R[aJ 1 \Q 2 ;f)} . 



In Fig.l, we give a concise summary of the IAA results for invariant analytic cou- 
plings a an (Q 2 , 3) and at(s, 3) calculated for one- , two- and three-loop cases in both the 
Euclidean and Minkowskian domains. 

Here, the dash-dotted curves represent the one-loop IAA approximations @ and 
(fPTl). The solid IAA curves are based on the exact two-loop solutions of RG equations^ 
and approximate three-loop solutions in the MS scheme. Their remarkable coincidence 
(within the 1-2 per cent limit) demonstrates reduced sensitivity of the IAA with respect 
to the higher-loops effects in the whole Euclidean and Minkowskian regions from IR to 
UV limits. 

For comparison, by the dotted line we also give a usual a s {Q 2 ) two-loop effective QCD 
coupling with a pole at Q 2 = A 2 . 

® 



As it has been shown in 21 



J, |25fl , relations parallel to eqs.(|15D and @ are valid 
for powers of the pQCD invariant coupling. This can be resumed in the form of a self- 
consistent scheme. Consider now new functional sets of nonpower perturbation expan- 



sions. 



6 As it has recently been established the exact solution to the two-loop RG differential equation for the 
invariant coupling can be expressed in terms of a special function W, the Lambert function, defined by 
the relation W(z)e w ^ — z with an infinite number of branches W n (z). For some details of analyticized 
and pipizated solutions expressed in terms of the Lambert function, see Rcfs. 31 [33| [l7| [34], [58). 
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Figure 1: Space-like and time-like invariant analytic couplings in a few GeV domain 



2.2 Expansion of observables over nonpower sets {.4} and {21} 

To realize the effect of transition from expansion over the "traditional" power set 

{a k s (Q 2 J)} = nj()-).ni. ...n 1 :... 

to expansions over nonpower sets in the space-like and time-like domains 

{A k (Q\f)} = a an (Q 2 J),A2(Q\f),A3... ; /)} = a(s, f), 2l 2 (s, /), 2t 3 . . . , 

it is instructive to learn properties of the latter. 
In a sense, both nonpower sets are similar 

— They consist of functions that are free of unphysical singularities. 

- First functions, the new effective couplings, Ai = a an and 2li = a are monotoni- 
cally decreasing. In the IR limit, they are finite and equal a an (0,3) = (5(0,3) ~ 1.4 with 
the same infinite derivatives. Both have the same leading term ~ 1/ In x in the UV limit. 

- All other functions ("effective coupling powers") of both the sets start from the 
zero IR values ^4fe>2(0, /) = 2tfc>2 (0, /) = and obey the UV behavior ~ l/(lna;) fc cor- 
responding to dig(x). They are no longer monotonous. The second functions A2 and 
2l 2 are positive with maximum around s,Q 2 ~ A 2 . Higher functions Ak>3 and 2tfc> 3 
oscillate in the region of low argument values and obey k — 2 zeroes. 
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Remarkably enough, the mechanism of liberation of unphysical singularities is quite 
different. While in the space-like domain it involves non-perturbative, power in Q 2 , struc- 
tures, in the time-like region, it is based only upon resummation of the "7r 2 terms". 
Figuratively, (non-perturbative !) analyticization in the Q 2 -channel can be treated as a 
quantitatively distorted reflection (under Q 2 — > s = —Q 2 ) of (perturbative) "pipization" 
in the s-channel. This effect of "distorting mirror" first discussed in [f27[] is illustrated in 
figures 1 and 2. 

Summarize the main results essential for data analysis. Instead of power perturbative 
series in the space-like D pt (Q 2 ) = 1 + d pt (Q 2 ) 



cWQ 2 ) = ]T4a s fe (Q 2 ;/) 

k>l 

and time-like regions Rpt(s) = 1 + r pt (s) 

r pt( s ) = 2^ r k a (s; /) ; (r 1)2 = di )2 , r 3 = d 3 - d x — - — , r 4 = d 4 - . . . ) , 

fe>i 

one has to use asymptotic expansions (|13"D and (|J) 

dUQ 2 ) = J2 dk MQ 2 , f) ; r v {s) = ^ 4 Ms, f) 

k>l fe>l 
with the same coefficients d^ over non-power sets of functions {A} and {21}. 



3 Global formulation of APT 

To apply the modified APT to analyze QCD processes, it is necessary to formulate it 
"globally", for the whole domain accessible to modern experiment, that is for regions 
with various flavour numbers / of active quarks. To this goal, one has to consider the 
issue of heavy quark threshold crossing. 



3.1 Threshold matching. 

In a real calculation, the procedure of the threshold matching is in use. One of the simplest 



is the matching condition in the massless MS scheme |39| 

a s (Q 2 = M 2 - f - 1) = a s (Q 2 = M 2 - f) (18) 

related to the mass squared M 2 of the f-th quark. 

This condition allows one to define a "global" function a s (Q 2 ) consisting of the smooth 
parts 



<W) = a s (Q 2 ; f) at M 2 < Q 2 < Mi 



(19) 
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and continuous in the whole space-like interval of positive Q 2 values with discontinuity 
of derivatives at the matching points. We call such functions the spline-continuous ones. 

At first sight, any massless matching, yielding the spline-type function, violates the 
analyticity in the Q 2 variable, thus disturbing the relation between the s- and Q 2 - 
channelsQ 

However, in the IAA, the original power perturbation series ( O ) with its unphysical 
singularities and possible threshold non-analyticity has no direct relation to data, being 
a sort of a "raw material" for defining spectral density. Meanwhile, the discontinuous 
density is not dangerous. Indeed, an expression of the form 

Pk(v) = Pk(o-, 3) + 6 (° ~ M f) {Phip\ f) - / - 1)} (20) 
with pk(cr; /) = Sa^-tx, /) defines, according to (|l4j) and flUf), the smooth global 



A k (Q 2 ) = - f Pk (a) (21) 

7T J a + x 
o 

and spline-continuous global 

oo 

SH*(a) = I ^-Pk(o-) (22) 

s 

functions 0. 

We see that in this construction the role of the input perturbative invariant coupling 
a s (Q 2 ) is twofold. It provides us not only with spectral density ( p0|) but with matching 
conditions (p~8|) relating A/ with A/ + i as well. 

Note that the matching condition (|18|) is tightly related [^, ^TJ to the renormalization 
procedure. Just for this profound reason we keep it untouched (compare with Ref. 



3.2 The s-channel: shift constants. 

As a practical result, we now observe that the "global" s-channel coupling a(s) and other 
functions 2lfc(s) generally differ from the effective coupling with a fixed flavor number / 
a(s; f) and 2lfc(s; /) by constants. For example, at Mf < s < Mf 

oo 00 

= y~p( a ) = J ' -^p( a '^) + y^-K 0- ; 6 ) =«(s;5) +c(5). 

s s Ml 

7 Any massless scheme is an approximation that can be controlled by the related mass-dependent 
scheme p0[ . Using such a scheme, one can devise [[ll) a smooth transition across the heavy quark 
threshold. Nevertheless, from the practical point of view, it is sufficient (besides the case of data lying 
in close vicinity to the threshold) to use the spline-type matching (|l^) and forget about the smooth 
threshold crossing. 

8 Here, by eqs.©,@ and @ we have introduced new "global" effective invariant couplings and 
higher expansion functions different from the previous ones with a fixed / value. 
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Generally, 

a(s) = a(s; f) + c(f) at M 2 < s < M 2 f+1 (23) 

with shift constants c(f) that can be calculated in terms of integrals over p(a; f+n) n > 1 
with additional reservation c(6) = related to the asymptotic freedom condition. 
More specifically, 

c(f - 1) = a(M|; /) - a(Mf] f - 1) + c(f) , c(6) = . 

These c(f) reflect the a(s) continuity at the matching points M 2 . 
Analogous shift constants 

Sl fc (s) = Ms; f) + c fc (/) at M| < s < M| +1 (24) 

are responsible for continuity of higher expansion functions. Meanwhile, C2{f) relates to 
discontinuities of the "main" spectral function (|20|) . 

The one-loop estimate with (3[^p{a\ f) = {ln 2 (a/Aj) + 7r 2 } , 

1 7T 1 7T X 7 f 

c(f-l)- c(f) = — — arctan - — arctan ^- ~ - a^(M|) (25) 

and traditional values of the scale parameter A 3 , A 4 ~ 350 — 250 MeV reveal that these 
constants 

c(5) ~ 3.10~ 4 , c(4) ~ 3.10" 3 , c(3) ~ 0.01 ; c 2 (f) ~ 3 a(M|) c (/) 

are essential at a few per cent level for a and at ca 10% level for 2l 2 • 

This means that the quantitative analysis of some s-channel events like, e.g., e + e~ 
annihilation f24 |, r-lepton decay p5| and charmonium width Q at the / = 3 region 
should be influenced by these constants. 



3.3 Global Euclidean functions. 

On the other hand, in the Euclidean, instead of the spline-type function a s , we have now 
continuous, analytic in the whole Q 2 > domain, invariant coupling defined, along with 
via the spectral integral 

oo 

"-W-Vj-^v** (26) 



with the discontinuous density p(a) (|20|). 

Unhappily, here, unlike the time-like region, there is no possibility of enjoying any 
more explicit expression for a an (Q 2 ) even in the one-loop case. Moreover, the Euclidean 
functions a an and Ak , being considered in a particular /-flavour region M 2 < Q 2 < 
Mj +1 , do depend on all A 3 , . . . , A 6 values simultaneously. 
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Figure 2: "Distorted mirror symmetry" for global expansion functions 

Nevertheless, the real difference from the / = 3 case, numerically, is not big at small 
Q 2 and in the "few GeV region" , for practical reasons, it could be of importance . 

This situation is illustrated by Fig. 2. Here, by thick solid curves with maxima around 
Q = A , we draw expansion functions A2 and 2t 2 m a few GeV region. Thin solid 
lines zeroes around A and negative values below, represent ^3 and 2l 3 . For comparison, 
we give also second and third powers of relevant analytic couplings a an and at . 

All these functions correspond to exact two-loop solutions expressed in terms of Lam- 
bert function []. 

4 Illustrations 

Another quantitative effect stems from the nonpower structure of the IAA perturbative 
expansion. It is also emphasized at the few GeV region. 

4.1 The s— channel 

To illustrate the qualitative difference between our global scheme and common practice 
of data analysis, we first consider the / = 3 region. 

9 Details of these calculations can be found in Ref.jlT]]. Assistance of D.S. Kurashev and B.A. Magradze 
in calculation of curves with Lambert functions is gratefully acknowledged. 
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The process of Inclusive e + e~ hadron annihilation provides us with an important piece 
of information on the QCD parameters. In the usual treatment, (see, e.g., Refs. 0, the 
basic relation can be presented in the form 

^ = l + r(s); r pt (s) = ^ + r 2 al(s)+r 3 al(s). (27) 

Here, the numerical coefficients ri = 1/tt = 0.318, r 2 = 0.142, r 3 = —0.413 (given for 
the / = 5 case) are not diminishing. However, a rather big negative r 3 value comes mainly 
from the —ri7r 2 /3? 5 J3 contribution equal to —0.456. Instead of fl27D , with due account of 
(H), we now have 

r„{s) = ^l + d 2 2l 2 (s) + d 3 2l 3 (s) ; (28) 

7T 

with reasonably decreasing coefficients d\ = 0.318; d 2 = 0.142; d 3 = 0.043, the men- 
tioned 7r 2 term of r 3 being "swallowed" by a(s)F] 

Now, the main difference between (|28|) and (^7]) is due to the term d 2 2l 2 standing 
in the place of d 2 a 2 . The difference can be estimated by adding into fl2"7| ) the structure 
r4« 4 with r4 ~ —1. This effect could be essential in the region of a(s) ~ 0.20 — 0.25. 
Here, in the APT analysis, the third, three- loop term contributes about half of a per cent, 
compared with 5,5% in the usual case. 

The APT algorithm with fixed / = 3 has recently been used [f43| for the analysis of 



Inclusive r-decay. Here, the theoretical expression for an observed quantity, the time-life 
of r lepton, contains QCD correction A expressed via an integral of an s-channel matrix 
element over the region < s < M 2 . 

As a result of the three-loop analysis of a modern [13J experimental value A exp (so 



3.16 2 ) = 0.191, it was obtained that a(M 2 ) = 0.380. Remind here that under usual 
treatment one obtains a s (M 2 ) = 0.334 that can hardly be related to any a s (M 2 ) value 
as far as the parameter 7r 2 /hi 2 (M^/A 2 ) is close to unity. 

Note also that the third term of (|28|) contributes here about one per cent. 



4.2 The Q 2 — channel : Sum Rules 

In the Euclidean channel, instead of power expansion like (0); we typically have 

d{Q 2 ) = + d 2 A 2 (Q 2 ) + d 3 A 3 (Q 2 ) . (29) 

7T 

Here, the modification is related to a non-perturbative power structures behaving like 
A 2 /Q 2 at Q 2 ^> A 2 . As it has been estimated above, these corrections could be essential 
in a few GeV region. 

10 This term contributes about 8.10 -4 to the r(M^) and, correspondingly, 0.0025 to the extracted 
a s (Mz) value. This means that the main part of the "traditional three-loop term" ^30:^ in the r.h.s. of 
(p^), being of the one-loop origin, is essential for the modern quantitative analysis of the data. In partic- 
ular, it should be taken into the account even in the so-called NLLA which is a common approximation 
for the analysis of events at y/s ~ Mz ■ For a more detailed numerical APT analysis of the / = 5 region, 
see g|. 
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In the paper f26|, the IAA has been applied to the Bjorken sum rules. Here, one has to 
deal with the Q 2 -channel at small transfer momentum squared Q 2 < 10 GeV 2 . 

Due to some controversy of experimental data, we give here only a part of the results 
of 261. For instance, using data of the SMC Collaboration |5| for Q 2 = 10 GeV 2 , the 



authors obtained a&x {Qq) = 0.301 instead of a^tiQo) = 0.275 . Here, the contribution of 
the third term is also suppressed. 

The same remark is valid in the analysis of the Gross-Llywellin- Smith (GLS) sum 
rules. Indeed, as it was shown in paper [p8fl , instead of proportions (65 : 24 : 11)tbE3 



of usual analysis, the APT gives (75 : 21 : 4)apt (for further details, see Section lie in 



|[28|| . The same effect for the Bjorken sum rules turns out[f2(J to be more pronounced 
(55:26: 19) TB ^ (80: 19: 1)atb- 

Some comments are in order: 

- We see that, generally, the extracted values of a an and of a are both slightly greater 
in a few GeV region than the relevant values of a s for the same experimental input. This 
corresponds to the above-mentioned non-power character of new asymptotic expansions 
with a suppressed higher-loop contribution. 

- At the same time, for equal values of a an (x*) = a(x*) = a s (a;*) , the analytic scale 
parameter A an values extracted from a an and a are a bit greater than that A^jg taken 
from a s . This feature is related to a "smoother" behavior of both the regular functions 
a an and a , as compared to the singular d s . 



4.3 Conclusion 

To summarize, we repeat once more our main points. 

1. We have formulated a self-consistent scheme for analyzing data in both the space- 
like and time-like regions. 

The fundamental equation connecting these regions is the dipole spectral relation (|l[) 
between renormalization-group invariant nonpower expansions D an (Q 2 ) and -R^-(s). 

Just this equation (equivalent to the Kallen-Lehmann representation), treated as a 
transformation, is responsible for non-perturbative terms in the Q 2 -channel involved into 
a an (Q 2 ) and non-power expansion functions {Ak(Q 2 )}- These terms, non-analytic in the 
coupling constant a, are a counterpart to the perfectly perturbative 7r 2 -terms effectively 
summed in the s-channel expressions a(s) and {21^ (s)}. 

2. As a by-product, we ascertain a new qualitative feature of the IAA, relating to its 
non-perturbativity in the Q 2 -domain. It can be considered as a minimal non-perturbativity 
or minimal non-analyticityQ in a as far as it corresponds to perturbativity in the s- 
channel. 

Physically, it implies that minimal non-perturbativity cannot be referred to any mech- 
anism producing effect in the s-channel. 

11 That is, the contribution of the first, linear in a s , is 65 % , while the contributions of the second 
and third are 24 and 11 per cent. 

12 Compatible with the RG invariance and the Q 2 analyticity — compare with p6| . 
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3. The next result relates o the correlation between regions with different values of 
the effective flavor number / . Dealing with the massless MS renormalization scheme, we 
argue that the usual perturbative QCD expansion provides our scheme only with step- 
discontinuous spectral density ( ptp depending simultaneously on different scale parameters 
Af ; / = 3, . . . , 6 connected by usual matching relations. 

This step-discontinuous spectral density yields, on the one hand, smooth analytic 
coupling a an (Q 2 ) and higher functions {Ak(Q 2 )} in the space- like region — eq. (|2T|) . 

On the other hand, it produces the spline-continuous invariant coupling a(s) and 
functions {2lfc(s)} in the time-like region — eq . (p2|) . 

As a result, the global expansion functions {Ak(Q 2 )} and {%lk(s)} differ both from 
the ones {Ak(Q 2 ] /)} and {2lfc(s; /)} with a fixed value of a flavour number. 

4. Thus, our global APT scheme uses the common invariant coupling a s (Q 2 ,f) and 
matching relations, only as an input. Practical calculation for an observable now involves 
expansions over the sets {Ak(Q 2 )} and {21a:(s)} , that is non-power series with usual 
numerical coefficients d& obtained by calculation of the relevant Feynman diagrams. 

In particular, this means that we have now three QCD effective couplings: a , a an - 
of the APT formalism, and traditional a s , . This usual one can be used for approximate 
expression of two first ones in four and five-flavor regions, for the comparison reasons. 

This means that, generally, one should check the accuracy of the bulk of extractions of 
the QCD parameters from diverse "low energy" experiments. Our preliminary estimate 
shows that such a revision could influence the rate of their correlation. 



5. Last but not least. As it has been mentioned in our recent publications [EI, 21 1, 
the IAA obeys immunity with respect to higher loop and renormalization scheme effects. 

Now, we have got an additional insight into this item related to observables and can 
state that the perturbation series for an observable in the IAA have better convergence 
properties (than in the usual RG-summed perturbation theory) in both the s- and Q 2 - 
channels. 
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